Complex WKB Evolution of Markovian Open Systems 
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We derive a semiclassical approximation for the evolution generated by the Lindblad equation as 
a generahzation of complex WKB theory. Linear coupling to the environment is assumed, but the 
Hamiltonian can be a general function of positions and momenta. The theory is carried out in the 
chord representation and describes the evolved quantum characteristic function, which gives direct 
access to the Wigner function and the position representation of the density operator by Fourier 
\ transforms. The propagation is shown to be of Liouville type in a complex double phase space, the 

. imaginary part of the action being responsible for decoherence. The theory is exact in the quadratic 

' case, just as the real WKB theory that we previously developed for the Markovian evolution of 

^N) I extended states, but it also describes the decoherent and dissipative evolution of localized states, 

such as the interference terms of a Schrodinger cat state. The present rederivation of the real WKB 
approximation leads to its interpretation as a first order classical perturbation of the complex theory 
and to a discussion of its validity. The example of a simple cubic Hamiltonian illustrates the various 
levels of approximation derived from the complex WKB theory. 



PACS numbers: 03.65.Sq,03.65.Yz 



Ph; I. INTRODUCTION 

^ \ The general state of a system in quantum mechanics can be fully described by the density operator p, also called 
^ ■ the "state operator", normaUzed by Trp = 1. If the system is known to be in a "pure state", that is, a single vector 
of the corresponding Hilbert space, then p is the projector on this vector, p = and then its Von Neumann 

I— entropy, S = — Trplogp, is equal to 0, and its purity, Tr'fP, is equal to 1. However, the density operator generalizes 
^ ' this case in that it can also describe a statistical mixture of such states, which expresses a partial knowledge of the 
, system. Then it is not a projector any more but a sum of several projectors, and one has S > and Tr'fp < 1. 
QQ • When the system is isolated, the evolution is unitary and preserves the purity of the state. However, when the 
QQ [ system interacts with its environment, which is always the case in realistic situations, an initially pure state will 
Cn| ■ undergo a non-unitary evolution which will not preserve purity: one is in the paradigm of quantum open systems. 

[ The consequent loss is two-fold. Decoherence corresponds to the vanishing of the off-diagonal terms of the density 
\^ • operator, leading finally to a statistical distribution which can be interpreted classically. Thus it can dynamically 
■ explain one aspect of the postulate of reduction of the wave packet. Dissipation corresponds to a loss or gain of energy. 
On , and is usually considered as much slower than decoherence, especially for large systems. 

' The general evolution of p{t) for an open system is a practical issue in many experimental situations, already well 
established in quantum opticsfl], also in atomic and nuclear physics, and more recently in the physics of quantum 
information where it is of crucial importance to controle the interaction of the system with its environment, in order 
to avoid decoherence. To construct theoretical models with the maximum of generality, it is generally admitted that 
' such an evolution should satisfy at least three basic requirements: preservation of Hermiticity {TrpA € TZ), trace 
(Trp = 1) and positivity {Trp\ip){ip\ > 0). 

Two main approaches have been explored to obtain explicit dynamics of an open quantum system. The first is to 
consider that the open system is part of a bigger closed system obeying a unitary evolution, and to trace over the 
environment part, that is over all the degrees of freedom which are not directly concerned with the system under 
study. Such an approach leads to the Bloch-Redfield equations The second is to assume that, beyond the 

above general requirements, the time evolution should obey a semigroup law, that is, a forward translational time 
invariance. Then the most general equation was derived by Lindblad[l| (see also 0]). 
The Lindblad equation can be written 



dp 
dt 



H,p 
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The Hamiltonian H describes the unitary evolution of the system without environment. There is no restriction 
on the operators Lk, commonly dubbed Lindblad operators, and they are not Hermitian when the coupling to the 
environment is dissipative. For instance, the master equation of quantum optics can be seen as a special case of the 
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Lindblad equation, with the operator Li = a, i.e. the annihilation operator, describing the emission process, and 
L2 = the absorption. 

In this paper, we will focus on finding semiclassical approximate solutions of a Lindblad equation with a general 
Hamiltonian part and Lindblad operators that are linear functions of momentum and position. This involves a 
considerable adaptation of the more familiar semiclassical approximations for pure states. In previous papers, we 
generalized the analysis of the evolution of extended states, the WKB theory based on a real classical phase space, 
that goes back to Van Vleck '6^ (see also [7|[8|f9(). Here we develop the WKB theory on a complexified phase space, 
thus generalizing previous semiclassical theory of unitary evolution developed by Huber, Heller and Littlejohn (lo| 
and by Maslov [nl . 

The advantage of the complex theory is that it allows us to include in the same theory density operators that 
are localized, such as, typically, the interference terms arising from superposition of coherent states. Also it is more 
robust to the fast damping of quantum correlations induced by decoherence, which breaks down the validity of the 
real stationary phase method. 

Just as with pure states, the evolution may in time extend the state beyond the range of validity for the theory 
and this intermediate stage also demands attention (see e.g. [3], for pure state evolution). However, this is a case 
where the treatment of open systems has an advantage: usually the process of decoherence is exponentially fast with 
respect to other relevant features of the motion. 

The complex WKB theory can also be used to rederive its real counterpart. Both become exact in the special 
case of a quadratic Hamiltonian and Lindblad operators that are linear in the positions and momenta. However, we 
find that the complex theory is more accurate, even for extended states for which the real WKB theory should hold. 
Fortunately, we are able to show that the perturbation theory that bridges both theories is farely accurate, if one uses 
the proper initial condition, which is fulfilled by a generic propagator introduced in [T^j- Therefore, it is not necessary 
to abandon the more transparent intuitive content of real WKB for the evolution of extended states. 

We represent the solution in the "chord" space, that is the Fourier conjugate of the Weyl-Wigner representation 
of the state operator, which is defined in the "centre" space. Together they describe a double phase space, where the 
double Hamiltonian generates corresponding classical trajectories that support the semiclassical approximation. It is 
this double phase space that is now complexified. The present treatement is in continuity of a series of papers, |il5|(l6j 
and [13], but it can be understood separately, since all the approximations are quite transparent. 

We first give a reminder for the Weyl-Wigner representation, and write the Lindblad equation in this formalism. 
Then we describe the general WKB procedure by introducing an adapted complex double phase space driven by a 
complex Hamiltonian, which leads to a Hamilton- Jacobi equation. After verifying agreement of the theory with the 
quadratic case, we develop an effective real WKB dynamics, based on a perturbative approach, which is valid for 
states whose initial classical action is real. We apply this real WKB to a generic mixed propagator Rx{y,t) which 
can turn an initial Wigner function Wo{x) into the evolving chord function xiy^'t)- Finally we construct a simple 
non-quadratic example where all the methods are compared. 

One should note that all the following formulae are appropriate for a system with a single degree of freedom, in 
order to clarify the notation. Nonetheless, it is quite simple to generalize our treatment and results for a finite number 
of degrees of freedom. 



II. LINDBLAD EQUATION IN THE WEYL FORMALISM 



The Weyl representation maps every quantum operator onto a phase space function, that is, a function of the vector 
X — {p,q)^Mi^M- -F'or an operator A the Weyl symbol A is defined as 

A{x) ^2jexp i-'-pQ) {q + - |> dQ. (2) 

The Weyl symbol of the state operator p is the Wigner function 

W{x)^Af J exp(--^pg) dQ, (3) 

with J\f = l/(27r?i); and its Fourier transform, the chord function x(^), also called characteristic function, is 

x{^)=Af Jexpi-^^Ax)Wix)dx, (4) 

where the wedge product of two vectors x = [p, q) and x' — {p', q') is defined by a; A a;' = pq' — p'q = Ix ■ x' , which 
also defines the skew matrix JJ = ( n ]■ One can have the chord symbol directly from the quantum operator 
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through the formula 



X(0 =AA J exp{-^^pq) {q+^\p\q- 
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dq. 



(5) 



We call the space of all x the centre space, and the space of all ^ the chord space. 

In the chord space, by using product rules for the product of operators [19|, the Lindblad equation is represented 
by a partial differential equation. This equation is actually simpler than in the Weyl (centre) representation, and this 
justifies our choice. In the case where the Lindblad operators are linear functions of p and q, that is, L = I' -x + il" ■ x 
with I and I real vectors, this equation can be written 



H{x' + l^,t) - H{x' - exp ( _ ^) A a;' ) x(l',t) d^' dx' 
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The dissipation coefficient, 



j = l"Al', 



(6) 



(7) 



is null for a Hermitian Lindblad operator {l" — 0) and we then have a purely diffusive case. H is the Weyl represen- 
tation of the Hamiltonian of the isolated system and coincides with the corresponding classical Hamiltonian, up to 
corrections coming from non-commutativity of p and q. Its arguments in equation ([6]) are the pair of remarkable points 

x+ = a; + S. and a;_ = x — ^, which can be considered as both tips of a chord ^. Although this chord ^ = (^p, can 
be interpreted as an auxiliary conjugate variable of x, in the current approach it is actually more convenient to write 
the solution in terms of y = = (— ^q,^p). Indeed, the direct sum of these conjugate spaces can be interpreted as 
a double phase space, where x formally plays the role of the position q, and y the role of its Fourier conjugate, the 
momentum p. Then the above equation becomes 



dx 
dt 



{y,t) = H{x' --Jy,t)-H{x' + -iy,t) exp [ -{y' - y) ■ x' ] xiv' ,t) dy' dx' 



{\"-yf x{y,t). 



(8) 



The same name has been kept for the characteristic function xiUi t), though strictly this should be x{^, t) — xi^^V^ t)i 
and we have set the complex vector 



X' + iX" = I{1' + il"). 



(9) 



The term y ■ ^ can actually be included in the integral term, thanks to an integration by parts of the exponential, 
and one has finally 



dx 
dt 



Af / U{x',y,t) exp -(y'-y)-a;' xt{y') dy' dx' 



1 



(A' . y)2 + (A" . y)2 xiv.t) 



(10) 



with 



'H{x,y,t) = H{x-^^y,t)- H{x + ^Iy,t)--f x-y 
= 'H+{x,y,t) -n^{x,y,t) - J X ■ y. 



(11) 



This is exactly the double Hamiltonian that generates the classical motion underlying the semiclassical approximations 
in p^l - Obviously, the double Hamiltonian will be time-independent if it is obtained from a time- independent single 
Hamiltonian. Furthermore, in the absence of dissipation, both Ti^{x, y) = H{x — -^Jy) and H~{x, y) ~ H{x + ^Jy) 
will also be constants which generate independent motions for both chord tips. 
Notice that (fTO|) can also be written as 

^[y,t)^~'-N I H{x',y,t) exp{-^yx')Wt{x') dx' 
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{X' ■ yf + {X" ■ yf x{y,t) 



(12) 



in terms of the evolving Wigner function. Alternatively, the definition of the complex double Hamiltonian 



Hc{x, y, t) = nix, y^t)~- [{X' ■ yf + (A" • yf 



(13) 



leads to 



(14) 



where and '•^^ mean that the derivatives are taken first and then the y terms are multiplied. Given this specific 
choice of ordering for the operators y = y and x = —j-§y, this is completely analogous to the Schrodinger equation, 
which thus allows us to extend the various forms of WKB theory, once it is recalled that variations in operator ordering 
have effects that are semiclassically small. 

The differential term in the RHS of (|14p (or the integral term in the RHS of represents the miitary part of the 
evolution. In other words, in the chord representation, the commutator is specified by the real double Hamiltonian as 



n 



chord 



h d 
i dy 



(1) 



x{y,t). 



(15) 



III. GENERAL COMPLEX DYNAMICS 



We here assume that the chord representation x(y, t) of the localized state has the usual semiclassical form 

x{y,t) = cxp-S{y,t), 



(16) 



where S{y, t) is a function with complex values of order 0{h^). One shall be aware that, if we find the time evolution 
of such a state determined by the Lindblad equation, then we can also evolve any linear combination of such states, 
which can be coherent states for instance. This is a consequence of the linearity of the Lindblad equation. 

This semiclassical form naturally induces a h expansion of the unitary part of the equation, as it is shown in the 
appendix \^ 



^\ i dy 



(1) 



y^^\t] x{y,t) 



n 



as, , 
-^{y.tly.t 



0{h) 



x{y,t)- 



(17) 



Hence, by expanding p4p at leading order in h, one obtains the Hamilton- Jacobi equation 



dS, , dS^ , 

-g^(y.t) = -n[-^iy,t),y,t 



-K 



iX'-y)' + iX"-yf 

dS_ 
dy 



0{h) 

{y,t),y,t] +0(h). 



(18) 



This leads to a double phase space generalization of the complex WKB theory in [l3|. Notice that the present 
Hamilton- Jacobi equation is defined by a Hamiltonian of unusual form: not only does it not separate into the familiar 
kinetic and potential energy terms, but it is complex. Therefore the action, S{y,t), becomes complex, even in the 
special case where the initial action, So{y) = S{y, 0) is real. Because of these unusual features, we will work directly 
with the Hamiltonian formalism in the full complex double phase space, rather than attempting to connect the 
evolving action to a Lagrangian. 

In general the initial action specifies a Lagrangian surface of half the dimension of the complex double phase space, 



yo 

Xo 



dSo . , 



(19) 
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To obtain the formal solution of the Hamilton-Jacobi equation, it is convenient to define a family of classical trajectories 
{xt,y^) in the complex double phase space with initial conditions on this initial surface and driven by the complex 
Hamiltonian l|13p through Hamilton's equations: 



Vt 



dn, ^ 



dx 

Xr = -^ixr,y^,t) = ^{xr,y^,t)-i{X' -y^) X' -i{X" -y^) X" 



(20) 



In the present case of a single degree of freedom, this family is spanned by variables z (2 complex dimensions) and 
T (1 real dimension). Hence, considering the real and imaginary parts of each complex variable, this forms a real 
submanifold M with 5 dimensions, within a real phase space of 8 dimensions, corresponding to the 4 complex dimen- 
sions of {x, y). This submanifold, which will serve as a backbone to build our solution, is completely parametrized by 
(y,r) through a function a;^(y,r), such that {xm{Vit)tV) is the most general point of M. 

In the real case, it is a well known fact, see for instance that the solution S{y,t) of the Hamilton-Jacobi 
equation 



^^{y.t) = -n[--iy,tu,t 



(21) 



is the generating function of the submanifold {xM{y,t),y) for each t. The complex generalization is straightforward 
and is proved in Appendix |B] by showing that the differential form 



Ss = -XMiy,T) ■ dy - Hc{xM{y,T),y,T)dT 



is a closed form on [y.T). Hence we can define the solution 



Siy,t)^Soiyi)+ Ss, 

Jy.^o 



y,t 



(22) 



(23) 



where y^ is arbitrary and the integral can be performed along any path of M. joining (y^, 0) to (y, t). 

A natural choice of such a path is (y^, t), such that (Xt, y^) is the classical trajectory of M with y^ — y, litterally 
the history of y at time t. We have in particular XMiyr^T) = by construction. The choice of this trajectory sets 
the value of y^ = yg- Then one can write more explicitely 



Siy,t)^So{yo) + 



t ^ 



dr 



lie {Xr,y^,T) 



dr. 



(24) 



One should keep in mind that (a;T-,y^) implicitly depends on {y,t), and in particular yg and Xq — —^^(j/o) ^re 
complex functions of {y,t). There may be many complex trajectories from the initial surface to that with fixed y, 
just as in ordinary complex WKB for unitary evolution [lol |. However, in the case that the initial action is real and 
the Lindblad terms are small, the relevant trajectory should lie close to a real trajectory of the unitary problem. 
This motivates a perturbative real WKB theory introduced in a further section. Tunneling might be included in the 
underlying classical description by allowing complex time trajectories. 

The final expression of the solution, in terms of the real part of the double Hamiltonian, is then 



I 1 
x(y, t) =K exp ( -Soiyg) ~ - 



_ 9y^ , _ V 

Xr ■ Vn{Xr,yr,T) 

OT 



2h 



(A'.y,)2 + (A".y,) 



dr 



(25) 



IV. THE QUADRATIC CASE 



Here we show that in the quadratic case, the solution ((25|) of the complex WKB approximation coincides with the 
exact solution derived in (Tsj . We start from a quadratic Hamiltonian 



H{x) = x-mx 

with some symmetric matrix H. The complex double phase space Hamiltonian, 



(26) 



Ticix, y)=-2x-my-jx-y-- (A' • y)^ + (A" • y 



(27) 
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induces the following double phase space dynamics: 

Xr = (2JH-7)a;^-i(A'-y^) A'-i(A"-y^) A" 

= (2EU + 7)y,. (28) 

Notice that in this case, even though is generally complex, this motion does not affect the trajectory y^, which 
remains real at all times. Then the "history" of (y, t) is 

Xr = -e-^^M.^(e-^*M7y) f [(A' • y JA' + (A" • y JA"] du 

Vr = C^^-'^Rl.y, ° (29) 

in the notation of section IIIII Here Rt defines the classical propagation operator generated by the Haniiltonian (|26p , 
that is 



exp(2JJH0, 



(30) 



and stands for its transpose. Because T-L{x,y,t) is linear in x, one has besides 

_ _ dy^ _ _ 

SO that, the expression of S{y^t) given by ([^^ boils down to 



(31) 



Siy,t)^Soie-"Rjy) + 



2 Jo 



dr. 



(32) 



Hence, the full complex WKB solution can be written as 



Xiy,t) = Xo(e '''Mf y) exp 



L / p27(t'-t) 



A • Rj_,,y 



dt' 



(33) 



Recalling our use of the same symbol, x(y), as a shorthand for x(^)^) = xi^^y^'t)^ together with we find that 
the complex WKB approximation coincides exactly with the exact solution in the quadratic case. Indeed, the general 
solution derived in [15i | was written in the following form 



Xit t) - Xo(e-^*M-t exp (^-^ e2^(*'-*) \l ■ Rt'-t$\ 



dt' 



(34) 



where xo($) is any general initial chord function. To see that ([M]) is identical with ([55]) one just need to notice that 



i-Jy 



(35) 



Thus, the general picture is that of exact classical propagation of the chord function (the same as the unitary 
evolution generated by the quadratic Haniiltonian), multiplied by a Gaussian factor that progressively attenuates the 
contribution of long chords. 



V. REAL DYNAMICS 



A. General analysis 



Complex WKB theory is more elegant than its real counterpart, in that the evolving complex action accounts 
explicitly for the entire evolution of the wave function, whereas the real WKB action has to be supplemented by an 
evolving amplitude. Even so, the advantage of working with a more intuitive real phase space may predominate, 
whenever both theories lead to similar results. It is important to note that, even though both sets of parameters A' 
and a" are real, the last integral in the general complex WKB solution ((25)) will not be real in the case of a complex 
trajectory. Thus it is not immediately clear that this still has the role of attenuating the amplitude of long chords 
as in the exact quadratic theory. In contrast, the real WKB theory which is also exact in the quadratic case, 
generalizes the quadratic exponent in (j34p by a negative decoherence functional with monotonical increasing modulus. 
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The necessity of the complex formahsm in the description of unitary evolution arises when the amplitudes, 0£j{^) 



(36) 



do not vary smoothly as compared to the complex oscillations. In this situation, simple (real) stationary phase 
approximations are not allowed, i.e., when the variation of real part of the action is also divided by the small 
parameter, h, it is necessary to resort to full complex saddle point approximations. Where this difficulty is not 
present, it was shown that both WKB theories lead to equivalent descriptions in the case of unitary evolution p^ . 
In short, the evolution of the imaginary part of the complex action matches that of the logarithm of the evolving real 
WKB amplitude. 

The formal similarity between the Schrodinger equation and the Lindblad equation (|14p for the chord function would 
permit us to immediately incorporate this equivalence, were it not for the imaginary part of the double Hamiltonian 
(fT3)) . Notice that the presence of the Lindblad dissipation coefficient, 7, in the real part of the double Hamiltonian 
(fTTj) is not a problem. The difficulty is that an initial real action cannot continue to be real, as it evolves classically 
due to a general complex Hamiltonian. The correct description of this feature is an advantage of the complex theory 
in the present context, which will now be examined. 

A simple rederivation of the real space theory developed in (l7j is to treat the imaginary part of the action through 
classical perturbation theory. In the absence of the imaginary term, a real action defines a real Lagrangian manifold 
through (Uni), which evolves through real trajectories. The first order effect of a perturbation of the Hamiltonian 
(whether real or complex) is to add a term. 



6Sp = —i 



Jo 



integrated over the (real) unperturbed trajectory, that we denote {xr^Vr) 
{xr,yT-). In the present case. 



(37) 

not to confuse it with the complex one 



Sn {x, y) = -i [(A' . + (A" • yf] = -1 |(A • ^/)|^ 



(38) 



so that exp{iSSp/h) is just the decoherence term in the real WKB theory [17]: It damps out the amplitude of long 
chords, thus cancelling the fine oscillations of the Wigner function. 

How crude is this approximation to the full complex theory? The way in which the classical perturbation is rederived 
in Appendix [Cl allows us to analyze this question. There, we consider a general one-parameter family of Hamiltonians, 
so that henceforth we replace SH{x,y) by a SH{x,y) such that, in the present case, the parameter in Appendix [Cl 
will be fixed as a = —i. Applying this exact general treatment of a one parameter family of Hamiltonians, the error 
in the perturbation theory, as compared to the exact action (IC4|) can thus be expressed as 



AP = 5S- 5S„ 



1 



da' / dr [5n{xr{a'),y^{a')) ~ 6n{xr{0),y,{0))] 



da' / dr-[i\'- yM')f - (A' • yA^)? + (A" • y.(a'))' - (A" • ^,(0))^] 
Jo ^ 



(39) 



But according to (PO)) . y is not directly affected by the parameter {a = —i), so that this error grows slowly as t is 
increased. Indeed, at any point in phase space the trajectory y^ has a velocity y that is independent of a, but the 
curvature of the trajectory is affected, because 



y{a) = m + aiX'-y)^X' + a (A" • y) A". 

Thus, keeping constant the initial variable jjoia) = y leads to the approximate trajectory, 

a<2 



y^{a) w y 



(A • y) -TTT A + (A -y) A 



dx"^ 



dx"^ 



and the initial estimate of the error is 

aH'^ 



AP: 



6 



[^■yf (A'-??A'U(A"-y)^ U'Z^X" 



dx"^ 



dx"^ 



(40) 



(41) 



(42) 
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This first correction of the complex action is entirely real, because the final parameter is = = — 1. If A = JHs 
considered small, as is usually assumed for Markovian evolution (the small coupling limit), then the error is OdZj*), 
while the decorerence functional is only 0(1^^). We also find that the decoherence functional grows linearly in time, 
while its error grows as t'^. A more subtle reason for the appropriateness of the real WKB approximation is that the 
decoherence functional damps the contribution of all large chords, so that, for long times, we are only interested in 
the classical region where y « 0. 



B. Real dynamics of the mixed propagator 

In and we introduced the mixed propagator Rx{y,t), such that 

x(y,t) = J Wo{x)Rx{y,t) dx. 
This propagator is a good candidate for applying the real dynamics because its initial expression is 



Rx{y,0) = exp 



(43) 



(44) 



which means that the initial action is purely real. Hence, the growth of the imaginary part of the action is initially 
perturbative, as will be verified in the following example. The important point is that this propagator can evolve any 
initial state, by using (^5]) . 

The previous perturbative approach gives the following expression for this propagator, 



2 2 

Rx{y,t) = exp ( --X -ilo + T 



dr-1- 

2h 



{\'-y^f + {^'-y^f 



dr 



where (xT,y^) is the real history of y at time t and with Xq = x, governed by Ti. and not Tic- Hence x^ and y^ are 
real functions of {x, y, t). 

This expression has the advantage of separating the Lindbladian damping term from the oscillating phase, which is 
not the case for a general complex WKB. Therefore one can conclude that for short enough times, the chord function 
is fading exponentially fast for values of y which are outside a disk centred in and with radius a/^, as it does in the 
quadratic case [l5!|. Moreover, one can also conclude from that it is also true for any initial state, although with 
possible modulation of velocity since the damping term depends on the initial x. 

In the next section we treat the case of an initially Gaussian state, as an example where the above calculated 
propagator can be used to get a general formula. It relies on the fact that the contributing x, in RxiUTt), ^^re then 
localized around the Gaussian's centre. 



VI. DYNAMICS OF A GAUSSIAN WAVE PACKET USING THE MIXED PROPAGATOR 

An initially localized wave packet will restrict the range of the required trajectories of (|25p. For instance, if one 
starts with a Gaussian wave packet centred in X, 

^o(a;) = ^exp(^-i(a;-X)'^, (46) 

then only trajectories whose starting point is within a distance ^/h from X will contribute to the integral. By 
expanding the generic trajectory around the central one {x^,y^) — {xt-{X ,y,t),y^{X ,y,t)), corresponding to the 
evolution of the the maximum X of the initial Gaussian, one can pcrfom the explicit integration of (|43p . 

If we call S{X,y,t) = S'{X ,y,t)+iS"{X ,y,t) the action such that i?j^(y,i) = exp jrS{X ,y,t), then Rj^^^,{y,t) 
can be expanded as 

By interpreting 5* as the generating function of the complex manifold driven by the complex Hamiltonian Tig, we 
show in appendix |D] that 

{X,y,t)^-yo, (48) 
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where i/q is the initial y coordinate of the complex trajectory (X,yQ,0) — > {xM{y,t),y,t), constrained to start at 
X = X and to stop at y in a time t. 
Then (gS]) gives 

„ = -1 / „p (-ix- + is(X, „, ,) - is. ^ + . fix') ,49) 

and then 

where is a complex valued fmiction of [X,y,t), hence a function of (y,t) once the initial state is fixed. 

The analysis holds as long as S{x,y,t) can be approximated by its quadratic expansion in x, that is, when the 
dynamics is sufficiently regular, or time sufficiently small. For a chaotic dynamics one expects exponentially diverging 
trajectories which may make expansion (I47p irrelevant for large times. 

If one uses the real WKB theory instead of the complex WKB theory, then the imaginary part of yp can be 
interpreted as a perturbation of the origin {/q of the real trajectory {x^ , y^): 

— (X,y,<)c.-y„ + * [X.y^)^\dT. (51) 
In the following example, both theories give the same result. 



VII. A CUBIC EXAMPLE 

Let us consider the special case of the Hamiltonian H{p, q) — and a linear Lindblad operator, L = Ip. The exact 
solution can be obtained, which allows us to check on our different levels of approximation. The Lindblad equation 
can be written 



dp i 
'dt 

In the p representation we then have 



2 



L{2ppp-fp-pp\) (52) 



|(p'ip(oip") = ~'-{p'''-p"'w\p{t)\p") - ^^{p'-p"f{p'\m\p"). (53) 



which gives the solution 



{p'\m\p") = {P'\m\p") exp ( -'-{p'' -p"')t-Ly -p")h]. (54) 



The general expression for the chord function is then 



X{y,t)^ I dp{p+y^\p{t)\p~y^)exp[--pyp 



= exp (^-^yl - t^yl^ J dp {p + ||p(0)|p - |) exp (^--^(y^ p + ty, 3/)) , (55) 

i.e. the original chord function is convoluted with an evolving Gaussian and then multiplied by another time dependent 
factor. 

The complex WKB approximation for this evolution is based on the classical dynamics ((20|) of the complex double 
phase space Hamiltonian (fTT|) . which is here equal to 



3 

ndp, q, yp, y,) = Sp% + ^ - '-{ly.f. (56) 
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It generates the following trajectories, 

iy_q)t = {yq)o = Vq 

{vp)t = {yp)o - GpoVqt 

3 

qt = qo + (3po + iJ/,)^ + ^l^Vqt 

Pt = Po- (57) 

One should note that in this simple example the chord motion, y^, is independent of the trajectory of the centre 
variable, Xt = {qt,Pt), just as occurs for quadratic systems. Hence, the evolved chord will remain real if it is real to 
start with. 

Given the initial action, 5o(y), its evolution is obtained explicitly from ([24]) as 

Siy,t)^So{yp + 6poyqt,yg)+t^3plyg-^+'-ily,f^ , (58) 
leading to the complex WKB approximation: 

i I ty^ it \ 

XwKBiy,t)=e^p- {So{yo)+tyq3pl-^ + -{lygf)\, (59) 

This can now be compared with the approximate saddle point evaluation of (|55p . Note that the initial density oper- 
ator in the p-representation, corresponding to (|16p with the initial action So{y), will also have a standard semiclassical 
form. It is just the Fourier inverse of (j55p for t = 0, with an action a^ij)' ^p") (a symmetrized Legendre transform of 
So{y)), such that 

^=-10 ^ d^=%- (60) 

Then the integral expression (|55p for the evolving chord function can be performed by the saddle point method. Thus 
the stationary point, po, of the exponent of (j55p is selected by the equation 



9o(Po + yg/2) - qoiPa - yq/2) = yp + Qty^pa = (2/p)o, (61) 

where, in the last equality, we have followed the yp-component of the classical trajectory in (j57p backwards in time. 
This trajectory depends on pq, a constant, according to (j57p . Hence, the full approximation for the evolving chord 
function becomes 



TTft / ityq tlyq i Vg 



Thus the phase of Xspiy^t) coincides with that of xwKB{y,t). 

In some cases this phase is exact. Consider, the simple example of a "Schrodinger cat state", composed of the 
superposition of a pair of coherent states, with the wave function 

1 f {q-Qf , ^Pq\ 



and different parameters, Xa — {Pa, Qa) and Xjj = (Pb, Qb). The density operator, corresponding to ji/'a) + \'4'b), has 
two diagonal terms and two off-diagonal terms: 

P= \\i>a)m + \\i^b)m + \\^a){M + \\i^b){i>a\. (64) 

We are interested here in the sum of the pair of non-local terms, Pab = IV'a)(V'6| and its complex conjugate. One can 
notice that the two first "classical" coherent terms can be retrieved by setting a — b. 
Inserting the p-representation of the initial density matrix, 

ip'lPabmP") = 7=^exp U{p'Qa-p"Qb) -^iP'- Paf - ^{p" - PbfX (65) 
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into (|55p . the integral can be performed to give the exact chord function evolution: 

_ 1 ( I AP^- {y^ + AQ)^ - AiPjy^, + AQ) 

-^(.,-APf-^^^%^^P^->-^%3), (66) 



where we have defined: 



P 



Pg+Pb 
2 

Qa + Qb 



AP - Pb-P, 

AQ = Qb-Qa- (67) 



The initial action obtained from setting i = is 



so that 



Soiy) = ~Qy, - ^(yp + AQ) + -{yp + AQY + -{y, - AP)', (68) 



dS i 
Po = --^(yo) = P+ -iyp + 6poyqt + AQ) 

90 = -|^(yo) = Q + ^(2/«-AP). (69) 



Then equation ([M)) gives 



S'(y, <) = S'o(yp + 6poy?i, y?) + t \ - ^ - ^('%)^ j 
- Pivp + Qpotyg) + |(yp + 6potyg + AQ)^ + 1 Uplvg - ^ - ^(^y,)' ) • (70) 



Here, 



1 + 32?;qt 

corresponds to the explicit evaluation of the stationary point defined by (|6ip . Therefore, both X5p(y,<) and 
XWKB{y,t) have the exact phase, being that the former also has the exact amplitude. 

This example of an initially localized wavepacket is a case where the real WKB theory becomes a crude approxi- 
mation, because the initial action is then complex and, hence, the contributing trajectory is also complex. 

However the simple classical nonquadratic evolution in this section does not lead to a complex trajectory in (j57p . 
if it is initially real, i.e., for an initial real action. In this case, the saddle point evaluation of (j55p becomes equivalent 
to the stationary phase method, leading to the real WKB result. In other words, the first order perturbation of the 
phase, which was studdied in section |Vl gives the full correction with respect to the complex WKB phase. 

In particular, the real WKB theory can be used to propagate the real initial action of the mixed propagator Rx{y, t) 
presented in section fVBl with its chord space representation (|44p . It corresponds to a density operator given by 

(P' + f Ipb' - f > = S{p' - p) exp [~'-q , (72) 
with X = {p, q) . Then, the integral in (j55p can again be performed exactly to obtain: 



it , I 



2 



III 



Rx{y,t) ^cxp[-—y^-t—yg^ exp[--yqq- -ypp- -tyqSp ]. (73) 



This expression coincides exactly with (j45p . which can be obtained explicitely from (j57p by setting 1 = 0. Thus, in 
this case both the real and the complex WKB evolutions are exact. 
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This is all the more interesting, in that Rx{y,t) can then be used to propagated any initial state. For instance, 
with the "Schrodinger cat state" , one gets 

Xab{y,t) = J Woix)Rx{y,t) dx, (74) 

with 

W^{x) = exp (^^^ + ^22/, , (75) 

which again gives This shows that the real WKB method, although inacurate for a general initial state, can still 
be used indirectly through Rx{y, t). 

Notice that expression ([74|) then also coincides with expression ((50|) . which is not so surprising as the action of 
Rx{y,t), in ((73|) . is actually quadratic in x. 



VIII. CONCLUSION 



We have shown that to leading order in ft, the chord function, or characteristic function, of a generic state will 
evolve according to a simple expression (j25p which is mainly determined by a generalized Hamiltonian dynamics 
taking place in a complex double phase space {x,y). The subspace y — corresponds to the ordinary classical phase 
space dynamics, whereas nonzero y accounts for quantum (nonclassical) correlations. 

When the classical trajectories used to build the phase of this expression are real, then the Lindbladian term is 
growing negatively and contributes to an exponential damping of the values of y which are beyond the disk of radius 
•\/ft around 0. This can be recognized as the region containing all the classical information about the state. The 
overall features of this evolution are certainly present in the case of a quadratic Hamiltonian and linear Lindblad 
operators, as derived for instance in [l3|, and fulfills our handwaving intuition that "decoherence drives the system 
back into classical dynamics" . 

We have here shown that this scenario also holds for general Hamiltonians, when the action of the initial chord 
function is real, and time short enough, so that the imaginary part of classical dynamics grows. Furthermore, it can 
be extended to any initial state by using the mixed propagator Rx{y,t), which fulfills the requirement for initially 
real action. 

For short times, the evolution of the initially localized components of a Shrodinger cat state, follow a simplified 
dynamics that may be understood without the full recourse to complex orbits. This is as true as of unitary evolution 
and the same jnethods employed in p^ . [l3l | can be generalized to a real double phase space. This approach has been 
presented in (23| 

For larger times the real and imaginary parts of the phase of the chord function become intertwined, and the 
evolution becomes difficult to unravel. However the decoherence time is known to be very small, and if complete 
decoherence is reached before the limit of perturbative validity, then the Wigner function becomes positive, and hence 
it will stay positive, because the Lindblad equation preserves this property. 



APPENDIX A: ASYMPTOTIC EXPANSION OF THE HAMILTONIAN OPERATOR 



One can verify easily that 



ld_ 

Xdp 



ld_ 
Xdq 



exp XS — 



dp 



dS 



d^S fdS 



dpdq \ dp 



1 / m{m- l)d'^S fdS 
^ Xy 2 

m — 1 / n o \ ^ 1 



m-2 



dS 



dS_ 

dq 



dp^ \ dp 

n{n-l) d^S fdS 



dq^ \ dp 



dS_ 
dq 



n-2\ 



-0{ 



A2- 



exp XS. 



From there one can generalize and write 
Ti' 



n 



dS_ 
dy 



,y 



2i ^ dx^ 



dS \ d^S. ,^,^2 



exp 1 -S 



(Al) 



(A2) 
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APPENDIX B: PROOF OF THE CLOSENESS OF THE DIFFERENTIAL FORM 5a 

We prove that the form 

5a{y,t) = -XM{y,t) ■ dy - n{xM{y,t),y,t)dt (Bl) 
is a closed form by showing that it obeys Schwartz equahties 

d d 

-Q^XM{y,t) = -^n{xM{y,t),y,t), (B2) 

and 

dpM{y,t) _ dqM{y,t) 



dyq dyp 

The point (xM{yTt),y) is on a particular trajectory (S^,y^) of the manifold A4 with 

xo = XM{y,t) 

iio = y 

(B4) 



(B3) 



and 



dn ,_ _ ^ , 

dy 
dl~L 

IJt = --g^ixT,yr^t + T) ■ (B5) 



One also has, of course, 



XM{yT,t + T) ^Xr- (B6) 

Now we can write the following sequence of equalities 

xm (y, t + dr)- xm (y, t) = XM{y,t + dr) - XM{y + ijodr, t + dr) + xm {y + y^dr, t + dr) - xm (y, t) 

= ^{y,t + dT)i/odT + Xr-xo + 0{dT^) 

dy 

dxM , , X dH ,_ _ , , dH ,_ _ , , , 9x 

^ ~Q:^^y'^^"-^''Q^^^o,yo,t)dT + -^(XQ,yQ,t)dT + 0(dT ) 

= ■^H{xM{y,t),y,t)dT + 0{dT^). (B7) 

This last equality shows (|B2p . 

On the other hand, at i = we have 

5p7w, „^ 9 dAo d dAo dqM , „^ 

oyq oyq dyp dyp dyq dyp 

dp 

which sets (|B3[) at i = 0. Now one just has to notice that -- — (y^, t + r) — (y^, t + t) is time independent: 

oVq dyp 

d (dpM,- , , s dqM,- , , A d d _ d d _ 
dr V dyq dyp J dyq dr dyp dr 

d dn,_ _ , d dn,_ _ 

= -J^-^iXr.yr.t + T) - — — {Xr,yr,t + T) 

dyp dyq dyq dyp 

= 0, (B9) 

which shows the statement. We used (IB5I) and (IB6I). 
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Therefore, since M is simply connected, 5a{y, r) is the exact differential of a function a(y, r) such that 



dy 

^''(y,T) = -n{xM{y.T),v,T). (BIO) 



d 

Hence one has 

so that A{y, t) = a{y, t) is the required solution. 

APPENDIX C: REAL WKB THEORY AS A PERTURBATION OF THE COMPLEX ONE 

Let us consider a familly of Hamiltonians, 7i(a;, y; a), where the parameter a may be real or complex. (It can also 
have more than one dimension.) It is then possible to reinterpret H as a single Hamiltonian in an expanded phase 
space: {x,y;a,f3) in which /3 is a hidden variable. It follows that a is a constant of the motion and the projection of 
each trajectory onto the {x, y) hyperplane coincides exactly with the trajectory of the original family of Hamiltonians 
for its specific value of a. 

Consider now the root search problem for finding, given any value of a, a trajectory that starts at the Lagrangian 
surface x — dSo/dy and arrives at specified value of y. Varying the parameter continuously in the interval [0, a], we 
obtain a one parameter family of trajectory segments, which form a finite 2-D strip in the enlarged phase space. This 
can be fully specified by choosing /? = as the arrival point for the extra coordinate (for all a) and its initial point 
to be 

f3o{a) = / dr $r ^ dr —{xria),y^{a);a). (CI) 
Jo Jo oa 

According to the Poincare-Cartan theorem [20| . 

j[x-dy + P-da-rLdT]^Q (C2) 

for any reducible circuit on a 2-D surface that is spanned by trajectories. So we can build this circuit from four 
segments: i. The trajectory that travels from the original Lagrangian surface and /3o(ck) to (y — const, (3 = 0) in the 
time t; ii. a path a 0, with all other extended phase space variables and the time held fixed; iii. the trajectory that 
reverses trajectory (i.), but with a = 0; iv. The circuit is closed (in zero time) by path along the original Lagrangian 
surface, while increasing a back to its original value, such that f3 = f3o{a), as defined by (jCip . Then (|C2p is rewritten 
as 

So{y{a))+ f [x{a)-ij{a)-H{Ar{a),y,{a))] dr 
Jo 

= Soiym + / [i(0) • ^(0) - n{Ar{0),yrm] dr - / /3o(a) da. (C3) 
Jo Jo 

Hence, the difference in the final action is exactly 

SS{a)^Siy,t;a)~S{y,t;0)^ P da' f ' dr ^{xria'),y,{a');a'). (C4) 

Jo Jo da' 

For small enough a, this may be approximated as 

AS{a)^ f dr [H(5;.(0),y,(0);a)-H(i.(0),y,(0);0)] =a / dr 5W(i,(0), y,(0)), (C5) 
Jo Jo 



if Hix, y; a) — Ti{x, y) + a STi.{x, y). This is just the first order classical perturbation: One integrates the change in 
the Hamiltonian along the unperturbed trajectory. 
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APPENDIX D: DERIVATIVE OF THE ACTION OF THE MIXED PROPAGATOR WITH RESPECT TO 

THE CENTRE SPACE VARIABLE 

For every Xq we define a manifold Ai{xo) by propagating the manifold a; = a;o in the double phase space driven by 
the complex Hamiltonian Tic- The coordinates of the manifold are (a;^(a;o, y, t), y, t), thus parametrized by {y,t). 
We want to evaluate the infinitesimal difference S{xq + dx^, y, t) — S{xo, y, t), where the variable Xq corresponds to 
the initial manifold, the variable y to the final observation point, and t to the time of evolution. We have 



S{xo,y,t) ^ -xo-yo- [ {xr ■ y^ -nc{xr,y^,T)) dr, 
Jo 



(Dl) 



where {xT,yr) = {x_\4{xQ,y^,T),y^) is the Hamiltonian trajectory which starts at (a;o,yo) time t = and ends 
at y at time t = t. On the other hand, one has 



S{xo + dxo,y,t) = -{xo + dxo) -yo- / XMixo,y,T)-dy + nc{xM{xo,y,T),y,T)dT, (D2) 

■^y<rO 



where we used the lagrangian nature of A4{xo + dxo) which gives us the freedom to start from the point (a;o + c?a;o, y^) 
at r = 0, which is indeed a point of the manifold, and follow any path which arrives at y on time t. In particular, one 
can chose the path {xm{xo + dxa,y^,T),y^,T), which follows "vertically" the A4{xa) trajectory, but which is not a 
Hamiltonian trajectory on A4{xo + dxo). Then one can see that 

S{xo + dxo,y,t) = S{xo,y,t) - dxo ■ y^ 

~J \^-g^{xo,yr,T) dxa-y^ + -^—{xo,y^,T) dxo- -^{Xr,yr,T)j dr 
= S{xo,y,t)-dxo-yQ, (D3) 



so that 

as 



^^{xQ,y,t)=-yQ. (D4) 
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